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We present the theory of bosonic systems with multiple condensates, unifying disparate models 
which are found in the literature, and discuss how degeneracies, interactions, and symmetries con- 
spire to give rise to this unusual behavior. We show that as degeneracies multiply, so do the types 
of fragmentation, eventually leading to strongly correlated states with no trace of condensation. 
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I. BOSE-EINSTEIN CONDENSATION, 
FRAGMENTED CONDENSATES, AND 
STRONGLY CORRELATED STATES 

Bose-Einstein condensation (EEC) is a very robust 
phenomenon. Because of Bose statistics, non-interacting 
bosons seek out the lowest single particle energy state and 
(below a critical temperature T^) condense into it, even 
though many almost degenerate states may be nearby. 
At temperature T — 0, the condensate contains all the 
particles in the system [l| . Although this remarkable phe- 
nomenon was originally predicted by Einstein for non- 
interacting systems [J, it became understood, starting 
with the work of F. London Q, that it occurs in strongly 
interacting systems such as '*He. While interactions re- 
move particles from the condensate into other states, the 
energy gain by macroscopically occupying the lowest en- 
ergy state (or another state, such as a vortex) is suffi- 
ciently great that the interactions in a Bose condensed 
system manage only to deplete a fraction of the conden- 
sate, but not destroy it This is the case for liquid 
^He as well as dilute gases of bosonic alkali atoms 6J . 

In certain situations, however ^a system does not con- 
dense into a single condensate 0, B 0. ^i ^> ^i 
Q, 0, 0, 0, 01 . In this paper we explore the physics 
of condensation when the ground state can contain sev- 
eral condensates - situations of fragmented condensation. 
This possibility arises naturally from the very concept of 
BEG ~ i.e., that the non- degenerate ground state of a sin- 
gle particle Hamiltonian is macroscopically occupied. 
How, instead, does condensation take place when the 
ground state of /lo is degenerate, with two or more states 
competing simultaneously for condensation? What hap- 
pens if the ground state degeneracy G is not just of or- 
der unity, but of order iV, the number of particles? And 
what happens if G becomes much greater than N , or 
even approaches infinity? How do the bosons distribute 
themselves in these competing levels? In all these cases, 
interaction effects play an important role in determining 
the structure of the many-body ground state. Different 
types of interactions produce different fluctuations (such 
as those of phase, number, or spin) and lead to different 
classes of ground states. 



Exact degeneracy is difficult to achieve, since small 
external fields or weak tunneling effects are generically 
present. Such energy-splitting effects lead to a non- 
degenerate ground state and therefore favor the forma- 
tion of a single condensate. Only when interaction effects 
dominate over energy splittings can a single condensate 
break up. 

Gonsiderations of the effect of ground state degenera- 
cies are not merely theoretical exercises. Rather, such 
degenracies occur in a wide range of current experiments 
in cold atoms. The cases where the degeneracy G is of 
order unity are related to bosons with internal degrees 
of freedom. Examples include a pseudospin-1/2 Bose gas 
made up, e.g., of two spin states \F = l,m = —1) and 
|F = 2, m = 1)) of ^■^Rb [l^, and a spin-1 or spin-2 Bose 
gas such as ^^Na or ^^Rb in an optical trap [2y|. In the 
former case, even though the two spin states of *''Rb are 
separated by a hyperfine splitting of order GHz, they can 
be brought to near degeneracy, with G — 2, by applying 
an external rf field. In the spin-1 Bose gas, the three spin 
states are degenerate at zero external magnetic field, and 
G = 3. The case of G ~ 7V^/^ is encountered in one di- 
mensional geometries where the density of states 
has a power law singularity at low energies. The case of 
G ~ is realized for bosons in optical lattices 22] with 
G sites each having a few bosons; in the limit of zero 
tunneling, G equivalent sites compete for bosons. The 
case G > > is realized in rotating Bose gases with very 
large angular momentum, L, in a transverse harmonic 
trap ^23;, !24j. As L increases, the rotation frequency 
of the atom cloud approaches the frequency of the trans- 
verse harmonic trap, wt, causing the single particle states 
to organize into Landau orbitals, which become infinitely 
degenerate as 11 — > lot- The great diversity of phenom- 
ena in these experiments is a manifestation of the physics 
of Bose condensation for varying degrees of degeneracy. 

As we shall see, when the degeneracy is low (G ~ 1), 
a single condensate can break up into G condensates. 
Following Nozieres 0, we refer to such a state as 
a fragmented condensate (define more precisely below). 
We caution, however, that there are many sorts of frag- 
mented states. Even degeneracy as small as G = 2 can 
give rise to fragmented condensates with distinct prop- 



2 



erties which depend on whether interaction effects favor 
phase or number fluctuations. System with larger de- 
generacies usually have a larger number of relevant in- 
teraction parameters and are more easily influenced by 
external fields. As a result, they have a greater vari- 
ety of fragmented states. For degeneracies G compara- 
ble to or much greater than the particle number, the 
system does not have enough particles to establish sep- 
arate condensates. Instead, interactions effects tend to 
distribute bosons among different degenerate single par- 
ticle states in a coherent way, establishing correlations 
between them. In such strongly correlated systems, inter- 
action effects obliterate all traces of a conventional singly 
condensed state. 

We stress at the outset that while interaction effects 
can cause fragmentation in the presence of degenerate 
single particle states, the presence of near degeneracies 
does not force the condensate to fragment. In many cases 
the ground state for a macroscopic Bose system is a con- 
ventional single condensate. Creating a fragmented state 
typically requires carefully tuning the parameters of the 
system, and whether such a state can occur in practice is 
very much dependent on the system. In optical lattices 
and double well systems, where the tunneling between 
wells can be tuned arbitrarily finely, a fragmented state 
can easily be achieved. Yet in other systems such as a 
spin-1 Bose gas in a single trap, or a rotating Bose gas, 
the parameter range allowing the existence of fragmented 
states scales like 1/iV, making it difficult to realize these 
ground states unless the number of particles is reduced 
to ~ 10"^ or fewer. 

In this paper, we focus on fragmentation in the ground 
state of Bose systems. We consider a number of canon- 
ical examples (double-well systems, spin-1 Bose gases, 
and rotating Bose gases) , which have increasing degener- 
acy in the single particle Hamiltonian. These examples 
illustrate the origin of fragmentation, the variety of frag- 
mented states, and their key properties. We will not 
discuss here fragmentation in optical lattices or in dy- 
namical processes, for they are a sufficiently large subject 
to require separate discussions. To begin, we first define 
condensate fragmentation more precisely, and the discuss 
general properties of certain classes of fragmented states. 



II. DEFINITIONS OF CONDENSATION AND 
FRAGMENTATION 

The concept of Bose-Einstein condensation was gener- 
alized to interacting systems by Penrose and OnsagerQ, 
|25| in the 1950s, by defining condensation in terms of the 
single particle density matrix, 

pW(r,r')-(^nO^(r)), (1) 

where V'^(r) creates a scalar boson at position r, and 
is the thermal average at temperature T . Since p'^' 
is a Hermitian matrix with indices r and r', it can be 



diagonalized as, 

p(i)(r,r') = 5]iV.(T)/,(r')*/.(r), (2) 

i 

where the Ni{T) are the eigenvalues and /i(r) the or- 
thonormal eigenfunctions of p^-^); / f*{r)fjr) — Sij. Set- 
ting r = r' and integrating over r we have Ni{T) — N, 
where N is the number of particles. We label the eigen- 
values in descending order, i.e., Nq > Ni > N2 > ■ ■ ■■ 
Equation implies that if one measures the number of 
bosons in the single particle state one finds Ni. This 
does not mean that the wavefunction of the many par- 
ticle interacting system is a product of such single par- 
ticle eigenstates. Unless there are special reasons (such 
as strict symmetry constraints, e.g., translational invari- 
ance), the eigenf unctions {fi} need not be the same as 
the single particle eigenstates {ui} of the single particle 
Hamiltonian of the non-interacting system. Two rele- 
vant examples are a homogeneous system of particles in 
free space, and a system in a harmonic trap. In the for- 
mer case, where the momentum k is a good quantum 
number, we have i — and {fi} and {ui} are the same 
plane- wave momentum eigenstates. In an inhomogeneous 
trapped system, there is no simple relation between the 
{fi} and the {uj} 

The usual situation of Bose condensation corresponds 
to the one eigenvalue A^o being of order TV, while other 
eigenvalues are of order unity, i.e., 

BEC: p(i)(r,r') - 7Vo(T)/o(r')*/o(r) + • • • (3) 

EE ¥°>{r')¥°\r) + --- , (4) 

or simply 

==7Vo|1'(°))(*(°)| + --- , (5) 

where 5''^°)(r) = -v/]Vo/o(r), and (• • •) denotes terms with 
eigenvalues Ni ~ 0{1). Since the macroscopic term in 
Eq. ^ is identical to the density matrix of the pure single 
particle quantum state, |^'^°^), the function \E'(°^(r) is of- 
ten referred to as the "macroscopic wavefunction" of the 
system. Systems in which p'-'^' has only one macroscopic 
eigenvalue, Eq. I^J, have single condensates. The ad- 
vantage of the Penrose-Onsager characterization of BEC, 
Eq. Q, is that it applies to both interacting and non- 
interacting systems, since it makes no reference to dy- 
namics. Penrose and Onsager also demonstrated the re- 
markable fact that Eq. Q holds for a Jastrow function, 
which is a reasonable approximation to the ground state 
of a system of hard core bosons, therefore substantiating 
Eq. I0J as a general property of interacting Bose systems. 

The Penrose-Onsager characterization can be easily 
generalized to bosons with internal degrees of freedom, 
labeled by an index /i. With field operator ipj^ir), the 
single particle density matrix is 

p«(r,M;r',p') ^ ir')¥^\r) + ■ ■ ■ . (6) 

In a conventional noncondensed system, such as a zero 
temperature gas of non-interacting fermions, or a high 
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temperature gas of bosons, all the occupation numbers 
are small: A^^ ~ ©(iV^). A conventional singly condensed 
system has one large eigenvalue A^o ^ C>{N), with all 
other eigenvalues small, A^i>o ~ 0(iV°). A fragmented 
system has q > I large eigenvalues, Ni^q ^ 0{N^). 
There is clearly a range of other possiblilites such as hav- 
ing an extremely large number q ^ 0{N^/'^) eigenvalues, 
each of which are of size N,^q - ©(iVi/S). This latter 
case occurs in an interacting system of one dimensional 
bosons, and is associated with a phase incoherent quasi- 
condensate [2lll27| . 



A. Simple example of fragmentation and its 
relation with coherent states 



in which TV bosons are condensed into the single particle 
state (e-^^/^a^ + e''*'/'^al)/V2. The coherent state is an 
example of a single condensate, where the single particle 
density matrix, p^^^ = {aj^a^^), is 

~ ^ \ e'*/2 ) ■ ^^^^ 

The difference between Eqs. (|13|l and H1U|I is the ab- 
sence of the off-diagonal term {a\a2) in the latter, which 
represents the loss of phase coherence in state (pseu- 
dospin) space. 



Before examining the origin of fragmentation in detail, 
let us consider a basic example of fragmentation - the 
Nozieres model j2| • Consider a system of N bosons each 
of which have available two internal states; 1 and 2. As 
we consider in detail later, this model can also be used to 
describe atoms in a double well potential. The Hamilto- 
nian of Nozieres' model consists solely of an interaction 
between atoms in the two spin states. 



H 



.9 t t 9 



(7) 



a]ai 



where the a] create bosons in state i = 1,2, and Ui 
is the number of particles in i. The interaction between 
the particles can be either repulsive [g > 0) or attrac- 
tive {g < 0). The eigenstates have a definite number of 
particles in each well |A^i, iV2), with N = Ni + N2, and 
energy 



E = -N1N2. 
2 



(8) 



Clearly, for g > 0, the ground state is two-fold degener- 
ate, with A^i = N,N2 = 0, or A^i = 0, iV2 = N; these 
states have single condensates, whose density matrices 
have eigenvalues and N. On the other hand, for g < 0, 
the state with Ni = N2 = N/2 has lowest energy. This 
i^ocfc-state, 



iN/2 JN/2 



(iV/2) 



|0), 



(9) 



has a fragmented condensate; the corresponding single 
particle density matrix. 



3 \a.f,a^) 2^0 1 



(10) 



has two macroscopic eigenvalues. 

One can contrast this fragmented state to the coherent 
state 
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1 



4^10) (11) 



B. Relation between Fock and coherent states 

The Fock state © is an average over all coherent phase 
states |0Ar), as we see from the relation. 



2N/2 J 



(14) 
(15) 



the latter relation holds for >> 1. As we discuss in 
the next section, this connection is very useful for un- 
derstanding the origin of various ground states. An im- 
portant implication of this relation is that for a macro- 
scopic system, the expectation value of any p-body oper- 
ator, Op ^ ojiifljis ■ ■ ■ '^lip'^i^p ■ ■ ■ '^1/2 '^I'l' ''^ Fock state 
is indistinguishable from that in an ensemble of coherent 
phase states \(t)), as long as p <C iV, 



(16) 



This equation, which we shall prove momentarily, 
shows that by measuring quantities associated with few 
body operators, one cannot distinguish a Fock state 
from an ensemble of coherent states with random phases 
|i2&i29j . An illustration of this effect is the interference of 
two condensates initially well separated from each other. 
Prior to any measurement process, the system is in a Fock 
state, Eq.lO, since there is no phase relation between the 
two condensates. Experimentally, in any single shot mea- 
surement (a photo of the interfering region), one finds an 
interference pattern consisting of parallel fringes whose 
location is specified by a phase 0, as if the two far away 
condensates actually had a well defined relative phase 
with this value [s^- The value of (j), however, varies ran- 
domly from shot to shot, so that if one averages over all 
the measurements, the interference fringes average out, 
as described simply by Eq. 31] . 

To prove Eq. (|16|1 . we use Eq. H15|) to write. 



ttN 



(17) 
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it is then sufficient to show that {(t)'j^\Op\(t>M) vanishes 
unless the phases and 4/ are very close to each other. 
We note that 

(0WI0iv) = cos^ ((0 - 0')/2) « e-^(^-^')V8 

^.J^^5{4>~4>'), (18) 

for N » 1. For p-body operators of the form Op = 



(^AflOpl^Ar) 

2P(iV-p)! 



g.[0(p/2-,)-^'(p/2-,')](0'^_^|^^_p).(19) 



For N >> we see from Eq. (HHJ that (0^|Op|0Ar) = 
v/8V^(0w|Op|0Ar)(5((/. - 0') + ©(l/Af), from which 
Eq. CSl) follows. 



tential [3^, internal hyperfine states coupled by electro- 
magnetic fields 20] , atoms in a rotating toroidal trap 
|33j . or wave packets in an optical lattice l34]. 

In solving this model it is useful to write Eq. (|20ll in the 
Wigner-Schwinger pseudospin representation |35j. We 
introduce the operators 

Jx — {a\a2 + alai)/2, Jy — {a\a2 — a\ai) /2i, 

Jz = - 4a2)/2, (22) 

which obey the angular momentum commutation relation 
[Ji, Jj] = isijkJk, i = x,y, z, and satisfy 



9 9 9 N /N 



(23) 



The Hamiltonian H2()|l can be written in terms of J as 

H = -2tJx + U{Jl +3"^ -N). (24) 



III. CHARACTERISTIC EXAMPLES OF 
FRAGMENTATION 



Mean-field solution 



We now develop three different examples that illustrate 
the origin of fragmentation. These examples are chosen 
to illustrate the increasingly complex behavior of frag- 
mentation when the number of degenerate single particle 
states increases. 



A. Scalar bosons in double well 

After the model of sec. Ill Al the simplest model with a 
fragmented ground state is that of bosons in a double well 
potential with tunneling between the wells. Unlike in the 
previous example, this model produces two distinct types 
of fragmented states. We label the wells by i = 1, 2; we 
assume that there is only one relevant state in each well, 
and that particles within a given well have an interaction, 
U , which can be either repulsive (U > 0) or attractive 
([/ < 0). We take the Hamiltonian to be 



H 



-t{a\a2 + 4ai) + ^ ["-i(?M - 1) + ?i2("-2 - 1)] , 



(20) 

where the aj creates a boson in well i, and Ui = a|ai is 
the number of particles in well i. The first term describes 
tunneling between the wells via a tunneling matrix ele- 
ment t (which we assume to be real and positive). The 
form Uni{ni — l)/2 — Ua\a\aiai/2 is the usual contact 
interaction {g/2) Jip^ip^ipip reduced to the single mode 
in each well. For fixed number of particles ni + n2 = iV, 
the interaction term can be written simply as 



(ni-na) +iV^-2iV 



(21) 



This model, simple as it is, has wide applicability to 
many physical situations: atoms in a double well po- 



As we shall see, the Hamiltonian (I2U|I can be solved 
exactly. Nonetheless, the mean-field solutions illustrate 
much of the physics of the true ground state. They also 
allow one to see the kind of fiuctuations about the mean- 
field state that lead to condensate fragmentation. 

The mean-field states are of the form of (pseudo)spinor 
condensates 



1 



{ua\ 



|0), 



(25) 



where u 



cos(6'/2) and v = e'''/^ sin(6l/2). The 
matrix elements of the density matrix in this state 
are {a\ai) = (ni) — N cos^{9 /2), (0^02) — (^2) = 
Nsm^le/2), and (ajaa) = iV sin(6l/2) cos(6l/2)e"^. The 
angles and (j) therefore characterize the density and 
phase difference between the bosons in the two wells. 
In the pseudospin language, the state H25|l describes 
a ferromagnet with total spin, (J) = (iV/2)n, where 
n — (sin cos 0, sin sin 0, cos 0) is the unit vector with 
polar angles {6,(j)). According to Eq. H24|l . its energy is 



= —tN cos 6 sin t 



(26) 



U ( — (C0S2( 



For a repulsive interaction, U > 0, E{6, (j)) is minimum 
at = 0, 6 — Ti/2\ orn = x. The mean-field approach 
therefore selects the non- interacting ground state |C) = 
{a\ + 4)^|0)/V2^iV! = |x) as optimal. Condensates in 
the neighborhood of |C) have energy 



i;(7r/2 + e, 0) = E{-k/2, 0) + ]^tN(t)'^ + ^ {2t+UN)e^ + . 



(27) 
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From this result one can begin to see problems with the 
mean-field solution: as t ^ with fixed U, the energy of 
phase fluctuations (0) vanishes; therefore quantum fluc- 
tuations begin to mix in many nearly degenerate phase 
states, \tt/2, 0). 

The mean-field solution for attractive U is very differ- 
ent from that of repulsive U. The solution depends on 
whether |;7|A^ < 2t or \U\N > 2t. In the former case, |C) 
is locally stable since the energy of fluctuations, described 
by Eq. H27|) . remain positive definite. However, the stiff- 
ness constant for density fluctuations (e) is lower than 
that for phase fluctuations (0). Thus, as \U\ increases, 
density fluctuations become dominant. Note, however, 
that the condition \U\N < 2t only occurs for a bounded 
number of particles, and under ordinary circumstances 
is not expected to be achievable in macroscopic systems. 
On the other hand, when \ U\N > 2t, we see from Eq. (|26f) 
that the optimal states satisfy sinfl = 2t/\U\N. There 
are two degenerate solutions: 9 = Oq = a.i-csm(2t /\U\N) 
and tt — Oq. For \UN\/2t >> 1, these two states approach 



we have 



- a}"|0)/VM and ~ 4"|0)/VM, corresponding to aU 
particles being in well 1 or 2 respectively. 



2. Exact ground states 

We now construct the exact ground state of the Hamil- 
tonian, (|20|l . and see how interactions can cause conden- 
sate fragmentation. At the same time, we can see how 
different types of interaction cause different types of frag- 
mented states. 

Non-interacting case: Let us first consider the simplest 
case of non-interacting bosons, with Hamiltonian H = 
—t{a\a2 + h.c). The single particle eigenstates are the 
symmetric state (ai -I- a2)/\/2 and antisymmetric state 
(ai — a2)/V2 with energy —t and t respectively. For a 
system of N bosons, the ground state is 



\C) 



1 



(4 + 4)^10) 



(28) 



with energy —tN. The single particle density matrix of 
this state is 




which has a single macroscopic eigenvalue X — N . The 
ground state is therefore a single condensate with con- 
densate wavefunction Cj = \/ N/2{1, 1) (the superscript 
"T" stands for transpose). 

Since the ground state |C) is a linear combination of 
number states |ni,n2) = a\^^ ai^^ |0)/ V^'i !?^2 ! , the num- 
ber of particles in each well fluctuates. 

We calculate the number fluctuations of the coherent 
state |C) by writing it in the number basis. For even N , 



\C) 



N/2 



(30) 



where \l) = |f +£,f -^), and 



(0) 



The number fluctuations arc then 



1/2 



-r/N 



(7riV/2)i/4 



(31) 



{AND = m {N,m 



d£ , = iV/4, 



(32) 

which, despite the approximations made in this deriva- 
tion, coincides with the exact result. 

Interacting Case: The many-body physics of this dou- 
ble well system is completely tractable. While one can 
calculate the properties the ground state numerically to 
arbitrary accuracy, we derive below all the essential fea- 
tures of the ground state by studying the effect of inter- 
actions on the non- interacting ground state, i.e., the co- 
herent state |C). We shall see that depending on whether 
the interactions are repulsive or attractive, the coherent 
state can be turned into one of two distinct fragmented 
states: a 'Fock-like' state or a 'Schrodinger cat-like' state. 

We first look at the Schrodinger equation of this sys- 
tem. Writing the ground state in the number basis, 



N/2 
f.=-N/2 



(33) 



we can write the Schrodinger equation, H\^) — E\^), 
where H is given by Eq. (|20|l . as 

E-^i = -ti+i-^i+i - ti^i-i + Uf-^i, (34) 

U = t^/{N/2 + £){N/2- e+1). (35) 

The many-body problem then reduces to a one- 
dimensional tight binding model in a harmonic potential. 
The special feature of this model is that the tunneling 
matrix element te is highly non- uniform [s^: te ~ A^/2 
for £ « 0, and te ~ \/N/2 for £ « ±N/2. This 
non-uniformity is a consequence of bosonic enhancement, 
a^N) = ^/N + l\N + 1) and a\N) = Vn\N - 1), which 
increases the matrix element by a factor of ^/N when re- 
moving a particle from a system with N bosons. As a 
result, {£ + l\a\a2\t) is maximum when both wells have 
equal number of bosons (i.e., £ ~ 0), and drops rapidly 
when the difference in boson numbers between the wells 
begins to increase, (i.e., £ » 1). A consequence is that 
hopping favors wavefunctions '^e having large amplitudes 
near £ = 0. For example, in the non-interacting case the 
wavefunction, Eq. (|31|l . is a sharply peaked Gaussian at 

e = o. 
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The interaction term, U{ni — 712)^/4, leads to a har- 
monic potential in Eq. Repulsive interactions {U > 
0) suppress number fluctuations, meaning that the Gaus- 
sian distribution [Eq. H31|l ] of the coherent state will be 
squeezed into an even narrower distribution. In the limit 
of zero number fluctuation. 



((Ani)^)-((An2)2)=0, 



the system becomes the Fock state 



1^) = ^ 



fN/2 fN/2 



{N/2) 



-|0), 



(36) 



(37) 



which is clearly fragmented, since it is made up of two 
independent condensates. This fragmentation shows up 
clearly in the single particle density matrix. 



/ 1 



= 2 V 1 



N f 1 
2 I 



(38) 



1 0) + ^(j )(0 1),(39) 



which has two macroscopic eigenvalues, corresponding to 
independent condensation in each well. The evolution 
from the coherent state |C) to the Fock state \F) can be 
captured by the family of states 



(7ra2/2)i/4- 



(40) 



As varies from 1 /N to much greater than unity, the 
initial coherent state |C) becomes more and more Fock- 
like. Indeed, exact numerical solution of Eq. (|34|) shows 
that Eq. (gnj is an accurate description of the evolution 
from a coherent state towards a Fock state. During the 
collapsing process when the wavefunction still extends 
over many number states (but fewer than 
one can take the continuum limit of Eq. H34|l . which re- 
duces to the equation for a particle in a harmonic oscilla- 
tor potential. From its Gaussian ground state wavefunc- 
tion, we extract a'"^ = {2/N)[l + UN/tf/"^. 

Using this continuum approximation, we calculate the 
off-diagonal matrix element. 



{aW) = 51 VW2~mW2T7TT)*^+i*£(41) 



« (iV/2)e-i/(2-'), 
which leads to a single particle density matrix 
N ( 1 



1 



The eigenvalues are Ai = -y(H-e ''^), and A2 
g-i/o"^) The relative number fluctuations are 

((Am)') = 0^12. 



(42) 

(43) 
-(1- 

(44) 



As 1/cr' varies from 1/iV to a number much larger 
than unity, the eigenvalues (Ai,A2) varies from (iV, 0) to 

varies from 



(iV/2,7V/2), and ^((AnT)') varies from ^J N to 0. 

This transition of the coherent state into a Fock state, 
described by the family Eq. (gUJ, is due to increasing 
phase fluctuations, as discussed in sec. nil A II The phase 
fluctuation effect can be seen by writing Eq. (|40|) in terms 
of phase states. 



E 



N 



|0) 



(45) 



where '^f'^ are the coefflcients given by Ea. (|31|l . The 
family of a. l40|) then becomes 



2x1/2 



(a7A^)i/4 



09 a^^^i 

2tt 



10) (46) 



(coherent state) to a^^ >> 1, the Gaussian in Eq. H46|l 
changes from a (5- function, 6{4>), to a uniform distribu- 
tion, thereby driving the coherent state towards a Fock 
state. 

Attractive Interaction: When U < 0, the potential en- 
ergy U{ni — 712)^/4 in Eq. (|21() favors a large number 
difference between the two wells, in particular, the states 
\£ = N/2) = \N, 0) and \i = -N/2) = |0, N). It therefore 
acts in the opposite direction as hopping, which favors a 
Gaussian distribution of number states around |^ = 0). 
The effect of interaction is then to split the Gaussian 
peak of the coherence state Eq. (|31|l into two peaks, a 
process which can be described by the family of states 



(47) 



where 2a is the separation between the peaks, cr' is the 
width of the peaks, and C is the normalization constant. 
As a varies from toward N/2, and a' shrinks at the 
same time from 1/Vn to 0, the state |^'(a)) = J2e ^'^l^) 
evolves from the coherent state |C) to a Schrodinger-cat 
state 



\Cat)^^{\N,0) + \0,N)). 



(48) 



The Schrodinger-cat state is fragmented in the sense that 
its single particle density matrix has two large eigenval- 
ues. 



{ala^) 



N 



1 
1 



(49) 



identical to that of Fock state. On the other hand, con- 
trary to the Fock state, it has huge number fluctuations. 



((Ani)2) = ((An2)2) = N^A. 



(50) 
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Details of how a and a' depend on the ratio U /t 
are found in |37j . where it is also shown that the fam- 
ily Eq. 147|) accurately represents the numerical solu- 
tion of Eq. (|34|l . This double- well example brings out 
the important point that fragmented condensates can- 
not be characterized by the single particle density matrix 
alone. Higher order correlation functions such as num- 
ber fluctuations are needed. This example also shows 
how a coherent state can be brought into a Fock state 
(or Schrodinger-cat state) through the phase (or number) 
fluctuations caused by repulsive (or attractive) interac- 
tion. 



B. Spin-1 Bose gas 

The spin-1 Bose gas, which is only marginally more 
complicated, provides an excellent illustration of the role 
of symmetry in condensate fragmentation. Here the de- 
generacies that give rise to fragmentation are due to a 
symmetry: rotational invariance in spin space. We will 
see that in the presence of local antiferromagnetic inter- 
actions all low energy singly condensed states break this 
symmetry. The true ground state, which is a quantum 
superposition of all members of this degenerate manifold, 
is a fragmented condensate. 

We consider a spin-1 Bose gas with 2N particles in the 
single mode approximation, where each spin state has the 
same spatial wavefunction. Several recent experiments 
have been carried out in this limit ^j .SS^ 4(i] (note that 
several of these experimental also explored regimes where 
the single mode approximation was not valid |4 ij). For 
conceptual clarity we consider the extreme case where the 
trap is so deep that all particles reside in the lowest har- 
monic state, trivially leading to the single mode limit. 
In this extreme situation there are no particle fluctua- 
tions: the only relevant interaction is scattering between 
different spin states. We denote the creation operator of 
a boson with spin projection = (1,0,-1) in the lowest 
harmonic state by aj^; the number of bosons with spin /i 
is Nfj^ = a)^a^. The most general form of the scattering 
Hamiltonian which is rotational invariant in spin space 




H = cS^ (51) 

where S — J2fj.u '^t.^fj-'^'^ij' with 5^^, the spin-1 matrices, 
{i = x,y,z), and c the interaction constant. We con- 
sider the case c > 0. More complicated Hamiltonians, 
with more interaction terms, are allowed for atoms with 
higher spins 43^ .4^ . The variety of fragmented states 
proliferate rapidly as the atomic spin increases. 

To illustrate how condensate fragmentation is affected 
by external perturbations, we include the linear Zeeman 
effect, with Hamiltonian, 

Hz = ~pS, = ~p{Ni - iV_i), (52) 

where p is the Zeeman energy proportional to the ex- 
ternal magnetic field B. To explain current experiments 



one also needs to include the quadratic Zeeman effect: 
the atomic energy levels of atoms are not linear in B due 
to hyperfine interaction between electron spins and nu- 
clear spins. We ignore these nonlinearities as they are 
irrelevant for describing fragmentation. 

For later discussion, we also include a term of the form 

Hq — e{a\a_i + a'^iQi), (53) 

which mixes spin states 1 and —1, where e is a con- 
stant. Terms of this form can be generated by magnetic 
field gradients in which case e is proportional to the 
square of the field gradient. Since both H and Hz con- 
serve Sz, the density matrix (aj^a^) for the ground state 
is diagonal in the presence of these terms alone, and the 
system is generally fragmented unless the density matrix 
happens to have only one macroscopic eigenvalue. The 
effect of He is to mix the 1 and -1 states, bringing the 
system into coherence, similar to the role of the tunneling 
term in the double well system. 

1. Mean-field approach 

Like the double well problem in the previous section, 
the Hamiltonian in Eq. (|51|l is exactly soluble. As previ- 
ously, we first consider the mean-field solution so that we 
can relate the fragmented condensate to a linear combi- 
nation of singly condensed states. 

The general form of a singly condensed spinor Bose 
condensate is 

Using the fact that a^a^lOjv = \/N{N ~ 1)C/3Ci'IC>a'^2' 
it is easy to see that the number fluctuations are {N'^) — 
(iV^)2 = (N,^) = iV|C^p. Writing the Hamiltonian as, 

H = S^^ • alal^apa^ + 2N, (55) 

we have, up to terms 0{N^^), 

{H + Hz)\o,^c{Sf -p{Sz) + 2cN (56) 

where (S) = ;y C^m" C- It is easy to see that 

Eq.(jnn|) is minimized by 

where, as before, T stands for transpose, <f> is the relative 
phase between the I and —1 component, iVi -I- A^_i = N , 
and the difference A/ = A^i — A^_i = {Sz) is given by 
the nearest integer to p/2c. Note that {Sx) — {Sy) = 0, 
and the entire family {C^} is degenerate, reflecting the 
invariance of Eq. (|56|l under spin rotation about the z- 
axis. In the absence of a magnetic field, p ^ 0, we have 
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Cl -> (e-^'^,0,e*'^) /y/2. However, since Eq. (ESJ be- 
comes fully rotationally invariant at zero field, any ar- 
bitrary rotation of the state (1,0,1) /\/2 is also an op- 
timal spinor condensate. The entire degenerate family 
(referred to as the "polar" family in literature) is given 
by 



(58) 




with energy 



polar 



2cN. 



(59) 



Before discussing the exact ground states, we introduce 
the "Cartesian" operators 

ai + a_i ai-a-i 
Ax = ^ — , Ay ^ — — , = ao. (60) 



V2 



iV2 



The important properties of these operators is that that 
under a spin rotation Uf^^a^, where U — exp{—i9 ■ 

S), A rotates like a Cartesian vector, i.e. Ai — > R{9)ijAj, 
where Rij{9) is a rotational matrix in Cartesian space 
(xyz). It is easily verified that [Aj, ^j] — Sij, and 



iA^ X A, iV = A^ • A, 



(61) 



(62) 



and the operator that creates a singlet pair is, 

0t = j^n = -2a\al^ + al^. (63) 
The mean-field state Eq. (j54|l can now be written as 

|a) = (a- At)^|0)/%/7V!, d*-d=l. (64) 

which has average spin 

(S) = -~iNa* X a (65) 

The optimal mean-field state is determined by minimiz- 
ing 

{d\H\a) = ciV^ [1 - |a • ap] - ipNz ■ a* x a. (66) 

In zero field (p — 0), the optimal mean-field state is 
a = fi, where n is a real unit vector. The polar family 
mentioned above can now be conveniently represented by 
all possible directions of n, and the polar state in Eq. 1)58(1 
corresponds to 



(n-At)^, , 
in) = ^ ^^|0). 



/TV! 



(67) 



2. Exact ground state in a uniform magnetic field: 

In the absence of a magnetic field, H — cS^. Since H 
is proportional to the angular momentum operator, the 
ground state for C > (with an even number of bosons 
N) is a singlet. The singlet state of a single mode spin-1 
Bose gas is unique '43^, and hence any manifestly rota- 
tionally invariant state must be the ground state. One 
can therefore obtain the ground state by forming sin- 
glet pairs, which for an even number of particles yields 



1^ = 0) oc (A^ • A^)^/2|Q^^ 



(68) 



Since this state is rotationally invariant, its single par- 
ticle density matrix (aj^a^) is proportional to the iden- 
tity matrix (Shur's theorem). With the constraint N = 
J2fi '^li'^fi^ then have 



(69) 



The condensate is therefore fragmented into three large 
pieces. The energy of the ground state is exactly zero, 
whereas that of the mean-field state is 2cN (as shown in 

Eq. isni)- 

To relate the exact ground to the optimal mean-field 
state (the polar family {|n)}), we note that 



IS* = 0) cx 



^(n-At) |0), 



(70) 



i.e., the exact ground state is an average over the family 
of optimal mean field states. This is the analog of the 
symmetry averaging relation in Eq. I|16|l . This averaging 
process represents the quantum fiuctuations within the 
family of degenerate mean-field states; the effect of these 
fiuctuations is to reduce the mean-field energy from 2cN 
to zero. 

Next we consider the case of non-zero magnetic field 
p ^ 0. the Hamiltonian is H + Hz = cS^ - pSz- The 
ground state is \S,Sz = 5") or simply denoted as l^*, 5), 
where S is an integer closest to p/2c. The state 15,5) 
can be easily obtained by changing a singlet pair into a 
triplet pair, and we have 



1 5, 5) =I?af (At2) 



.{N-S)/2 



|0), 



(71) 



where 2? is a normalization constant, and Bose statistics 
require S to be even if N is even. Since the Hamiltonian 
H + Hz conserves 5z, the single particle density matrix 
remains diagonal, {ajj^a^) = N^5^v As shown in ,1L.1 j|. 
the eigenvalues Nfj_ have a very interesting behavior as a 
function of S, namely 



TVi 



N{S + l) + S{S + 2) 
25 3 



(72) 



N-.^^-^^S^. No^^. (73) 



25 -f 3 



25 + 3 



9 



If S is of order 1, all three eigenvalues, N±i,Nq, are 
macroscopic. On the other hand, if S becomes macro- 
scopic, (i.e., S/N is less than but of order 1), both A^i 
and N^i remain macroscopic (with iVi — > {N + S)/2 
and A^_i — > {N — S)/2) while A^o becomes of order unity. 
This means as S increases, the Sz — component is com- 
pletely depleted. Although the system is still fragmented, 
the number of fragmented pieces is reduced from three 
to two, even for a tiny spin polarization. Further anal- 
ysis fl3| shows that as one increases S, the fluctuations, 
{an}) = (Nf) ~ {Ni^, drop rapidly. 

The reduction in the number of fragmented pieces mir- 
rors the previously discussed reduction in the size of the 
space of degenerate mean- field states. In the absence of a 
magnetic field, polar states with arbitrary n are degener- 
ate, while a magnetic field in the z direction favors those 
with n pointing in the x-y plane. 

To connect the exact polarized states to the mean-field 



states, we use the relation A^^ 
to write 



-2alal^- 



^0 ' 



Eq. 



(At2)(^-s)/2|o) 



(N-S)/2 
p=0 



S-2p,p) (74) 



where the exact form of Dp is not important for our 

discussion. If we now act on (A''''^)(^~'^)/^|0) by a\, 
the coefhcients of the large p states in Eq. H74|l will 
be bosonically enhanced via a||A^i) = {Ni + l)\Ni + 
1). Further acting on the same state by a\ eventu- 
ally picks out the term with largest p, which is | {N + 
5)/2,0, {N-S)/2), in the sum Eq. {HJ. This process is 
equivalent to replacing 4'^(-2a|ali -f aj)^)(^~'^)/2|0) by 
al^(-2alali)(^"^)/2|0). We then have 

\S,S) oc af(At2)(^-5)/2|0)«af(-2alali)(^-^'/2 



t(Af+S)/2 UN-S)/2 



|0). 



(75) 



From our previous analysis of the two-well system it is 
clear that this state is formed from an angular average 
of Eq. (|57|l . The field gradient term 1)53(1 plays the role 
of tunneling in the two-well system and can drive the 
system from a fragmented to coherent state. 

To complete the connection between these examples, 
we note that when the occupation of the Sz = state is 
negligible, we can set oq = 0, and only two spin states are 
required to describe the system. Under these conditions, 
the Hamiltonian for the spin 1 gas [Eqs. (|^ . 
reduces to the two- well Hamiltonian H2UI) with U — AC, 
and t = — e, and an asymmetry between the wells given 
by p. This analogy breaks down when p and the 
Sz — state becomes occupied. 



C. Fast Rotating Bose Gas: 

In the previous cases of pseudospin-1/2 and spin-1 
Bose gases, the number of degenerate single particle 



states is of order unity, far fewer than the number of 
bosons. Bose gases with large amounts of angular mo- 
mentum have just the opposite behavior; the number 
of nearly degenerate states can be much larger than 
the number of particles. As we shall see, for relatively 
small angular momentum, fragmentation similar to that 
in pseudospin-1/2 and spin-1 Bose gases can occur. How- 
ever, for very large angular momentum, it is energetically 
more favorable for a repulsive Bose gas to organize itself 
into a quantum Hall state, which is an extreme form of 
fragmentation in which all traces of conventional conden- 
sation disappear. Before discussing the fragmentations of 
rotating Bose gas, we first discuss how the single parti- 
cle energy levels of a rotating Bose gas turn into Landau 
levels and achieve high degeneracy as the angular mo- 
mentum of the system increases. For simplicity, we limit 
the discussion to two dimensional systems. 



1. Lowest Landau levels (LLL) and the general properties 
of many body wavefunctions in LLL: 

The single particle Hamiltonian for a particle in the 
rotating harmonic trap is 

ho-^Lz = + T:MLjlr^ - r!z • r x p. (76) 

where uj± is the trap frequency, and is the rotational 
frequency of the trap. Equation (|7()|l can be written as 



— — + -M(uj^. -n^)r^. (77) 

2M 2 ^ J \ J 



If one rewrites p — Muj±z x r as p — eB x r/2c, with 
eB/Mc — 2lu±, one sees that Eq. H77|l is identical to 
the Hamiltonian of an electron in a magnetic field Bz in 
a harmonic potential with reduced frequency uj'j_ — il^. 
To diagonalize the single-particle Hamiltonian, Eq. I|7t)|l . 
we begin by noting that the 2D simple harmonic oscil- 
lator Hq = p2/2M -I- MLu]_r^/2 is diagonahzed as ho — 
hu!±{al.a^ + alytty + 1), where — (d/h)px — ix/d,ay = 
{d/h)py—iy / d, d^ = h/mLu±, and Lz = —ifi{al.ay—aj^a^). 
Defining a± = {ax ± iay)/^/2, we have 

Lz — — ?i(a^a_|_ — a^_a_), (78) 

and 

ho—^lLz = h{LJ±+n)a''_^a_^_+h{i^J±—n)a''_a_+huJ±. (79) 
The eigenstates are therefore 



a\ aL 



(80) 



with eigenvalues 



En,m = ?i(w± + + h{LL!± — Q)m + huj±, (81) 
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FIG. 1: (Color online) Energy eigenvalues of a two dimen- 
sional harmonic oscillator as a function of angular velocity 
n. For clarity only seven values of m are shown, so that the 
infinite degeneracies at f2 = appear to be only seven-fold. 
The inset shows the energy states at f2 = as a function of 
their angular momentum L. (Colors: each Landau level is 
colored with a distinct hue.) 

where n,m — 0,1,2,.... The energy levels, shown in 
Fig. ^ are organized into Landau levels, labeled by n, 
separated by SEi = h{LL!j_ + Q). States within a given 
Landau level are labeled by m with spacing dE2 = h{uj± — 
f2). At criticality, Q = a;_L, each Landau level becomes 
infinitely degenerate. 

The many-body Hamiltonian, including a contact in- 
teraction g, is 

N 

H = Y,{h,-^lL,^)+gY,5{v,-Vi) (82) 

i—l i>j 

As ^ uj±, mixing is strongest among the states in the 
same Landau level. To simplify matters, we consider only 
extremely weak interactions, gn « 2hfl. This limit is 
naturally reached when is close to uj±, where the cen- 
trifugal potential largely cancels the trapping potential 
and the cloud becomes large and dilute. For such weak 
interactions, only the lowest Landau level is populated. 
The eigenfunctions Mn=o,m(r) in the lowest Landau level, 

^„,(r) = ^^e-NV2<i^^(r|^), z ^ x + zy, (83) 
V Tra^m! 

are angular momentum eigenstates with L ^ m, and a 
spatial peak at ^ {m\r'^\m) = d\j2m + 1. 

The many-body wavefunction for a systems of N 
bosons is then 

Umi(ri)...u™„(rjv), (84) 

[m] 

where {m} denotes the set of non-negative integers 
(toi, m2, . . . m^r). Since, apart from the Gaussian fac- 
tor, the single particle wavefunctions are of the form 



the many-body wave function is 

*(ri,...,rjv) = /(zi,...,zjv)e-S.=il^»lV2'i^ (85) 

where / is an analytic function which is symmetric in z^. 
In particular, a Bose-condensed state corresponds to 

JV 

^'u,(ri, ...,rjv) = ]^w(ri), w(r) = y^amttm(r). 

i=l m 

(86) 

If is also an eigenstate of the total angular momentum 
Lz = L, then / in Eq. H85|l must be a homogeneous sym- 
metric polynominal of {zi} of degree L. Moreover, the 
single particle density matrix in the angular momentum 
basis must be diagonal, i.e., 

(aLa«) = ^m«(4iam>- (87) 

Thus a Bose condensed state (Eq. (jHU) that is also an 
angular momentum eigenstate must be fully condensed 
into the single particles state u™, i.e., w(y) = Um(i"). 
The total angular momentum must then be L = mNh. 
Such a coherent state, however, offers no flexibility to re- 
distribute particles to lower the energy, and is therefore 
unlikely to be the ground state at fixed angular momen- 
tum or fixed rotational frequency f2, except for certain 
special cases in which parameters of the system are care- 
fully tuned. A general Bose condensed state, Eq. H8t)|l . 
would be one where w{r) contains more than a single an- 
gular momentum state, Um- Such states break rotational 
symmetry and are not angular momentum eigenstates. 
Furthermore, due to rotational symmetry, any rotation of 
w (u;(r, (/)) — > w{r, (j) -\- 9)) will have the same energy and 
angular momentum, L = {^w\Lz\^w)- One can there- 
fore form an angular momentum eigenstate with angular 
momentum L by forming the average 

$L(ri,..rjv) = y"d^e*^^*({r„(/., + 0}). (88) 

Such averaging always lowers the energy of the system, 
and generically results in a fragmented state. Specific ex- 
amples follow. Note that with small modification these 
arguments also apply in the general case when interac- 
tions are not weak. 



2. Ground state of the rotating Bose gas with attractive 
interactions: 

For the rotating Bose gas with attractive interaction, 
Wilkin, Gunn, and Smith 9] pointed out the ground state 
with non-zero angular momentum is one with all angular 
momentum carried by the center of mass, and has the 
form 

vl'(ri,..rjv) = if(Z/d)^e-I^"il^'l'/2'^', (89) 
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where Z is the center of mass and is a normaUzation 
constant, 



in which 



1 ^ 



K = 



(90) 



\N m! 



(99) 



It is straightforward to show that the single particle den- 
sity matrix of this state is 



(aL«ri> = ^rnnN 



L\{N - 1) 
m!(L - my.N^ 



From this degenerate set of coherent states we construct 
an eigenstate of angular momentum by taking the super- 
position 



(91) 
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The ground state is fragmented because Eq. has 
a distribution of large eigenvalues. Later, Pethick and 
Pitaevskii pointed out that this state is of the form 



J 27r 



/L!(^rd2pv 



(100) 



vl/(ri, ..r^) = ^cuimreiiim (92) ^here ^ = ^ E^i ^^/d 



where vI/cm(R) = {Z/d)^e^^^^''/'^''^ is the wavefunction 
of the center of mass R — J2iLi ^i/^^ (f>rei{{pi}) — 

ife~^i=i'^'' ^^"^ is a product of single particle states, 
where the particle coordinates are {pi — Vi — R}. Given 
this structure, it is natural to refer to this state as being 
singly condensed in the center of mass frame. Here, we 
show that this exact ground state can also be written 
as a symmetry average of broken symmetry states of the 
form Eq. H88|l . as in many of the previous examples. 

Attractive interactions favor particles clumping to- 
gether. In homogeneous systems, such clumping leads 
to collapse. In the fast rotating limit, however, the an- 
alyticity of the wavefunction in the lowest Landau level 
and the Gaussian factor impose strong constraints on the 
degree of localization possible. The most localized state 
is uo(r), a Gaussian of a width given by the trap length d. 
A similar localized wave-packet at location a = -\- idy 
is 



which gives 



(93) 
/V^d, (94) 



(95) 



A many-body coherent state formed from these single 
particle states is 



N 



(96) 



This state carries angular momentum hL = ('I'|Lz|vI>) 
= (*| /dr-0l'(r)r X (W/i)-0(r)|*), with 



L = Nj dr\Ur)\' {[r/d]' - l) = N\a/d\' 

The single particle density matrix is 

(Vnr')^(r)) =^0:(r')0a(r), 



(97) 



(98) 



3. Fragmentation in a rotating Base gas with repulsive 
interactions: 

Similar fragmentation is found in the repulsive case, 
which for L <C Nh is described by a vortex lattice. A 
general mean field state is of the form Eq. H86|l. with 
w(r) = /(z)e-l^l'/2'i', 



f{z) = Y[{z~a), 



(101) 



where the zeroes a are the locations of the vortices p7l |. 
The angular momentum carried by this state can be ob- 
tained by noting that within the lowest Landau level, 



{L^)^h{{ir/df-l)), 



(102) 



which increases as the density of vortices increases. Fol- 
lowing the arguments which we used for the attractive 
case, we are once again lead to a fragmented condensate 
(see ^3 for numerical studies). 

The rotating Bose gas with repulsive interactions is, 
however, much richer than that with attractive interac- 
tions, since the mean-field picture breaks down in a fun- 
damental way at large values of angular momentum. One 
sees this breakdown by first noting that with increase of 
the angular momentum, the vortex density increases and 
the particle density decreases 0,'49,'5^. Eventually, the 
density of vortices is comparable to the density of parti- 
cles and one can significantly improve the energy of the 
system by correlating the positions of the vortices with 
the positions of the particles. These correlations cannot 
be captured by any simple manipulation of the mean- 
field states; the general description of the system with 
large angular momentum is quite complicated ^3, 24] . 

Despite this complexity, there is a limit in which we can 
find the exact ground state. Imagine that the system is 
rotating with Q sufficiently close to lu± that the energy 
spacing within the lowest Landau level may be treated as 
a perturbation. One would find the ground state by first 
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minimizing the interaction energy in the lowest Landau 
level, and then perturbatively including the level spacing. 
For short ranged interactions the energy is minimized by 
any state for which the wavefunction vanishes whenever 
two particles come together. Degenerate perturbation 
theory then says that the lowest energy state of all these 
wavefunctions is the one with lowest angular momentum. 
In the lowest Landau level, the lowest angular momentum 
bosonic wavefunction which vanishes when two particles 
touch is the = 1/2 Laughlin state, 

^i[r])^Y[iz,~z,fe~\^^\'/^^\ (103) 

i>j 

For an infinite system, one readily sees that 

{alaj = (l/2),5„„ (104) 

for all TO. Not only are there no eigenvalues of order 
N, but they are all less than unity. All traces of the 
conventional Bose condensation are obliterated. 



IV. DISCUSSION 

A. Salient features 

The examples shown in Sec. Ill share a number of com- 
mon features. Many of the fragmentation processes we 
have discussed, such as those in Eqs. lO, and 

(|100|l . can be described by a family of quantum states 
that are weighted averages of broken symmetry states 
over the space of broken symmetry, typically of the form 

|*(A)) = J dxW{x,\)\x)c , (105) 

where \x)c is a coherent state with broken symmetry pa- 
rameter X (6-g-, a spin direction), W is sl distribution 
function in the space of broken symmetry, and A is the 
parameter that controls the fragmentation of the system. 
If, as A changes, say, from to 1, changes from a dis- 
tribution sharply peaked at xo to a completely uniform 
distribution, Eq. (|105|l will evolve from the coherent state 
Xo to a fragmented state. Such changes in the distri- 
bution function reflect a growing fluctuation about the 
initial coherent state xo singled out by a tiny symmetry- 
breaking field in the coherent regime. 

For example, in the case of repulsive Bose gas in a dou- 
ble well, X is the relative phase between the condensates 
in the two wells and xo = 0. In the spin-1 Bose gas, x is 
the vector n and xo is a direction normal to an infinites- 
imal external magnetic field. In the rotating attractive 
Bose gas, x is the location in space of the coherent state, 
and Xo is the equilibrium location determined by, e.g., a 
weak potential that breaks rotational symmetry. Exactly 
how the fluctuations about xo grow depends on the spe- 
cific dynamics of the system. Large fluctuations about 
Xo directly reflect the competition different degenerate 



states for Bose condensation, the ultimate cause of frag- 
mented structure, as discussed in the Introduction. 

Another important feature of states of the type H105|) 
is that, as a consequence of bosonic enhancement, the 
range of the control parameter A over which the system 
switches from a coherent to a fragmented state shrinks 
with particle number N. For a repulsive Bose gas in 
a double well, as discussed in Sec. III. A. 2, the transi- 
tion from a fragmented to a coherent takes place around 
t/U > l/N [33 • A similar situation occurs in the spin- 
1 Bose gas, where the singlet (fragmented) state give s 
way to a coherent state for field gradient G > l/N [l3| . 
Since the window for fragmentation vanishes as TV in- 
creases, fragmented condensates will not be realized in 
condensed matter systems, with TV ~ 10^^ particles. On 
the other hand, in mesocopic systems like quantum gases 
with typically N ~ 10^ particles, fragmented states can 
exist in parameter ranges accessible to experiment. 

Although our discussions have mainly focussed on the 
class of states Eq. (|105|) . other fragmentation processes 
do not lead to states that are most naturally expressed in 
this manner, e.g., the Schrodinger-cat family in Eq. H47|l . 
One can have distinct fragmented states, such as the Fock 
state and the Schrodinger-cat state, with entirely differ- 
ent properties but identical single particle density matri- 
ces. Such different states are therefore indistinguishable 
within the Penrose-Onsager scheme. To tell them apart, 
it is necessary to examine second order correlation func- 
tions such as number fluctuations. Even higher order 
correlation functions are needed to distinguish the frag- 
mented states of more complicated systems, such as those 
of F = 2 Bose gases of ®^Rb, and the recently realized 
spin-3 Bose gas of ^^Cr ^J. In general, one can expect 
a large variety of fragmented states, differing from each 
other by high order correlation functions. 

The relative number fluctuations of a fragmented state 
(such as those in the pseudo-spin 1/2 and spin-1 Bose 
gases discussed in Sec. Ill) are a measure of the sta- 
bility of the state. Huge fluctuations, such as in the 
Schrodinger-cat state and the singlet state of the spin- 
1 Bose gas indicate that the system is easily damaged by 
external perturbations. Consider, for example, a pertur- 
bation of the form H' = rjc^ai -l-h.c, where ai is a boson 
in one of the wells in the double well example, or a boson 
in the Sz = 1 spin state of a spin-1 Bose gas, ct adds 
a particle in a different atomic state of the same boson 
(e.g., a plane wave state in the background gas), and 77 
is a very small parameter. The perturbation H' acting 
on the Fock state a\^^'^al^^^\0) simply changes it to an- 
other Fock state, c^a|^^^ ^Oj^^^jO), which also has zero 
number fluctuations, (iV^ — {Ni}"^) = 0. In contrast, H' 
acting on the Schrodinger-cat state (a|^ + a2^)|0) col- 
lapses it into the state c'''a|'^~"^|0), immediately reducing 
the enormous number fluctuations of the Schrodinger-cat 
state to zero. More generally, if the cat state l^'cat) = 
1^1) + 1^2) is a sum of two Gaussians in number space, 
as shown in Eq. (|47|) . with one Gaussian (l^*!)) peaked 
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at |A''_|_,iV_) and the other 



1*2)) peaked at |A''-,iv + 



, the action of 
and suppresses 



where N+ + N_ = N, and N+ » N_ 
H' on l^'cat) considerably enhances 
1*2). 

The large effect of the small perturbation H' is due 
to bosonic enhancement, which gives H'\N^,N_) = 
^JN+ ct|iV+ - l,NJ) and H'\N_,N+) = JnZ ct|iV_ - 
l,N+). Since N+ » iV_, the norm of H'\N_,N+) 
is considerably smaller than that of H'\N+,N^). As 
a result, il'|*cat) « -ff'l^i) ~ x/^ct|iV+ - l,iV_), 
which is no longer a Schrodinger-cat state. In a spin- 
1 Bose g similar action changes the singlet state 
to one with much smaller relative 

number fluctuations. 

Since Schrodinger cat-like states can collapse into Fock 
states with the slightest perturbation, and Fock states 
can easily be reassembled into a single condensate by 
any small amount of tunneling between different frag- 
mented pieces, why should one bother with fragmented 
states? Is fragmentation relevant? The point, as men- 
tioned before, is that even though fragmented states can- 
not be realized in macroscopic systems, the situation is 
different for mesoscopic systems like trapped quantum 
gases. The huge reduction in particle number consider- 
ably relaxes the constraint of formation of Fock states 
and Schrodinger cat-like states, and fragmented ground 
states become realizable. The phenomena of fragmenta- 
tion becomes even richer if it takes place in both real and 
spin space, such as with high spin bosons in an optical 
lattices close to the Mott limit. The combined effect of 
spin degeneracy and spatial degeneracy (due to differ- 
ent isolated wells) produces a great variety of quantum 
phases as the spin of boson increases. In addition, the sin- 
glet ground state of a spin-1 Bose gas can also be viewed 
as a "resource" for singlets and may therefore be useful 
in developing quantum teleportation protocols in optical 
lattices with spin-1 bosons. Our understanding of the 
properties of fragmented states, and in particular their 
dynamics, is at such an early stage that it leaves open 
considerable room for inventive ideas, which is where the 
excitement lies. 

We stress that we have described only the simplest 
types of fragmented states. Many other kinds which do 
not fall into the category we considered. The v = \/2 
Laughlin state discussed in Sec. III.C.3 is an example at 
the other extreme. There are no large eigenvalues singled 
out in the density matrix. Instead, all the eigenvalues 
are identical and of order unity. In addition to level de- 
generacy, temperature effects can cause condensates to 
fragment, a subject to which we turn in the next section. 



B. Role of temperature 

Up to this point, the examples of fragmentation given 
in this paper have focussed on the case where interaction- 
driven quantum fluctuations break up the condensate. 
It is intuitively clear that thermal fluctuations can play 



a similar role. A trivial example is given by a non- 
interacting Josephson junction, governed by the Hamil- 
tonian 



H = -t{a^b + h^a). 



(106) 



This Hamiltonian is diagonal in the basis of symmetric 
and asymmetric states, for which the creation operators 
are (a^ -|-6^)/\/2, and (a^ — 6^)/-\/2, so the single particle 
density matrix is diagonal in this basis. At zero tempera- 
ture only the symmetric state is occupied - a single con- 
densate. At very large temperature each of these states 
is equally occupied so the system is fragmented. For in- 
termediate temperatures, the occupation of each of these 
states is 



2tl3m 



p2tl3m 
l^m=~N/2 ^' 



(107) 



N 



coth(/3t(iV + 1)) - ^ coth(/3t) 



The upper signs denote the symmetric state (s) and the 
lower the antisymmetric state (a) . The crossover between 
the singly condensed state at T = and the fragmented 
state at T ^ t is smooth. Extensions of this argument 
are relevant for spinor condensates in which one can, in 
principle, have a hierarchy of transition temperatures, 
where the fc = mode becomes macroscopically occupied 
below some temperature Tc, but order in t he sp in-channel 
does not occur until a lower temperature |53| . 

A vortex lattice provides a qualitatively different ex- 
ample of how finite temperature fragments a condensate. 
Imagine a bucket of ^He rotating at frequency fl. The 
ground state of the system contains a triangular array of 
vortices with vortices per unit area. At finite temper- 
ature the vortex lattice is thermally excited, giving rise to 
a decay in the phase correlations across the sample. If we 
let ipi''') be the superfluid order parameter coarse-grained 
on a scale l arg e compared to the vortex spacing, then ac- 
cording to ^3113, the correlation function {tp* {r')'tp{r')) 
decays as |r — r'\~^, for large separations. The expo- 
nent f] = l/(37r^/cirvA^) is proportional to the ratio of the 
distance between particles to the distance between vor- 
tices. Here the particle number density is p, the thermal 
wavelength is A = lmk\,T and the distance between 
vortices is r^. This algebraic decay of correlations in real 
space corresponds to an algebraic decay in momentum 
space. For ry < 3, the occupation of the fc = mode 
scale as A^o ~ L^~^, and as as fc ^ the occupations of 
the fc 7^ modes scale as Nk ~ k^~^. Thus the number 
of macroscopically occupied modes scales as L'', and for 
carefully chosen L and rj a fragmented state can result. 

For a typical helium experiment 77 < lO"'^, and the 
depletion caused by this effect is negligible. Experiments 
on alkali gases, e.g., create vortex lattices with 

Tv ~ 5/im and a particle density p ~ 10^^ cm~^, for which 
r] ~ 10"'^, also a minor correction. Experiments with 
smaller vortex lattices in Paris 1561 have a similar vortex 



14 



spacing ~ 2/xm, and density p ^ 10^'^ cm"^, yielding 
a comparable value for rj. Although vortex lattices in 
current experiments are not thermally fragmented, there 
does not appear to be any fundamental impediment to 
making rj larger. 



C. Other recent work 

Not surprisingly, given the fascinating nature of frag- 
mented states, several recent papers (all cited in the 
relevant sections of this paper) have been investigating 
models similar to those discussed here. As a guide to 
the reader, we briefly summarize these works, presenting 
them in roughly the same order as they appear in the 
main text. 

The two-state system was explored by Nozieres 
who pointed out the existence of fragmentation in that 
system. The generic stability of a two-state fragmented 
condensate was considered by Rokhsar [Tsl l . One realiza- 
tion of this model would be to place atoms in bonafide 
double-well potentials. Due to mode mixing such a 
system is not identical to a simple two-mode model, 
and several authors have explored more realistic mod- 
els. Spekkens and Sipe J6] used a variational approach 
to compare fragmented and singly condensed states in 
such potentials. More recently, Streltsov and Cederbaum 
[s?! . along with Moiseyev ,58] produced similar results 
through a multi-mode mean-field theory. The evolution 
of the ground state from a Fock state to a coherent state 
and from a coherent state to a Schrodinger-cat state was 
studied by Ho and Ciobanu [s^l- Experiments have re- 
centl y p roduced a condensate confined in a double-well 
trap |33 |. 

The fragmented nature of the ground state of the spin- 
1 Bose gas was first noted by Nozieres and Saint James 
0. Later studies by Law, Pu, and Bigelow 45] investi- 
gated the ground state properties and the spin dynamics 
in terms of the basis used in Eq. (|74|l . Ho and Yip [l^ 
showed that the fragmented singlet state has huge fluctu- 
ations; they also showed the relation between this frag- 
mentation and spontaneous symmetry breaking. They 
explicitly showed how magnetic fields and field gradients 
drive the system into a singly condensed state. Similar 
considerations were addressed by Koashi and Ueda [l]| . 
Javanainen discussed Issues involving the measurement 
of the fragmented spin-1 ground state [T^ . 

The rotating attractive gas was first studied by Wilkin 
et al. who noted that the ground state was frag- 
mented. The connection between this fragmentation and 
symmetry breaking was first discussed by Pethick and 
Pitaevskii 0. 

There have been extensive studies of the properties of 
the rotating repulsive gas in the lowest Landau level - 
many of which have focussed on the structure of the or- 
der in the single particle density matrix. These works 
have mainly used a combination of exact diagonalization 
and variational techniques. A particularly relevant pa- 



per is the exact diagonalization study of Liu et al. [43 , 
which focusses on the symmetry-breaking nature of the 
vortex states. A related paper by Jackson et al. [S^ 
compares mean-field and exact wavefunctions at various 
values of the angular momentum. The Goldstone mode 
associated with vortex nucleation was discussed by Ueda 
and Nakajima (6C|. 

Fragmentation occurs in "clumped" bosonic systems 
with attractive interactions. Ueda and Leggett used a 
two-mode approximation |33j | to study fragmentation and 
soliton formation in a one dimensional attractive Bose 
gas. A quite thorough comparison of mean-field the- 
ory and exact diagonalzation is found in the articles by 
Kanamoto et al. [61| , as well as the closely related work of 
Kavoulakis 0| . Montina and Arecchi [ssj use a Monte- 
Carlo scheme to investigate the degree of fragmentation 
of this system. In three dimensions, Elgaroy and Pethick 
[T^ showed that a harmonically trapped gas of atoms 
with attractive interactions does not form a stable frag- 
mented state. 

Other systems with symmetry breaking and fragmen- 
tation include phase separated two-component gases 64| , 
and rotating gases during a phase-slip event (3^. Bo- 
son ground states where the condensate is broken into a 
macroscopic number of pieces include the Mott insulator 
"251 , fractional quantum Hall states psL , and the low 
dimensional Bose gas (27^ ,65j . Aspects of the dynamics 
of regaining phase coherence among condensates are dis- 
cussed by Yi and Duan 0|, a paper closely related to 
discussion of how the measurement process is influenced 
by fragmentation [28lf29|. 

Finally we mention the lecture notes produced by 
Castin and Herzog for the 2000 Cargese Summer School 
[Tsl , which lucidly introduces fragmentation and analyzes 
the case of spin-1 bosons and of the one-dimensional at- 
tractive gas. They included an extended discussion of 
the role of symmetry breaking. 



V. SUMMARY 

Bose condensation is remarkably robust. It is therefore 
exciting to search for zero-temperature bosonic states 
which are not condensed and to understand the process 
that breaks up the condensate. Fragmentation, where 
the condensate breaks up into a few pieces, is the first 
step in this journey, which eventually ends at strongly- 
correlated states possessing no trace of condensation. 

We have described some of the canonical models which 
have fragmented ground states, and extracted their prop- 
erties. We see that there is a rich variety of fragmented 
ground states. Often fragmentation is associated with 
restoring a broken symmetry. Sometimes it is accompa- 
nied by order in higher-order correlation functions. The 
one unifying feature appears to be a combination of near- 
degeneracies and interactions. The higher the degree of 
degeneracy, the more fragmented the condensate may be- 
come. 
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